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Abstract 

(^ We examine in pure SU{3) the dependence of extended monopole 

^—^ current k and cross-species extended monopole current fc'^™'^'^ on tem- 

(^ perature r, monopole size L, and fractional monopole charge oc 1/q. 

"^ We find that features of both k and k^^°^^ are sensitive to r for a range 

5^ of L and q. For example, the spatial-temporal asymmetry ratios of 

j^ both k and k^^°^^ are sensitive over a range of L and q to the SU(3) de- 
confinement transition. The motivation for studying cross, extended. 



^H and fractionally charged monopoles in SU(3) is explained. 



1 APQCD 

A challenging problem in QCD is to identify its confinement mechanism and 
understand how it works. To this end "compact" or lattice QED(CQED), 
whose action is 

- ScQED = Yl I^CQED COS 0p, (f ) 

P 

provides a compelling prototype. Upon a BKT transformation [f], the CQED 

expectation value of a Wilson loop W = exp z(A, J) in lattice different forms 

notation is 

(W/)oc J2 exp{-^}, (2) 

{k\dk=o} 

S = —^ (J, A-V) + 27r' l3cQEDik, A-'k) - 2^i{dk, A-'*E). (3) 

2pcQED 

J and k are, respectively, conserved electric and magnetic monopole current 
loops. E is the electromagnetic field due to external current J: dE = J. The 
first and second terms of S correspond to the electromagnetic interaction 
energies of J and k. The third term is the interaction between the monopole 
currents and the background electric field E created by J. When coupling 
constant ficQED is sufficiently small, the entropy of the sum over monopole 
loops in (2) dominate over suppression by Boltzmann factor exp{ — S*} and 
monopoles condense. In this phase CQED exhibits the dual Meissner effect: 
monopoles dynamically rearrange the background electric field E into an 
effective fiux tube [2] which gives the Wilson loop its string tension. In this 
way, magnetic monopole condensation causes electric confinement in CQED. 
An analogous demonstration that monopole condensation is the origin 
of QCD confinement would be a great achievement. To this end 't Hooft 



suggested the following idea [3]. Suppose QCD monopoles, like the 't Hooft- 
Polyakov monopoles of the Georgi-Glashow model [4], carry charges that are 
magnetic with respect to the [t^(l)]^~^ Cartan subgroup of color SU{N). 
Then SU{N) gauge symmetry obscures the magnetic charges and it is nec- 
essary to gauge fix at least the SU{N)/[U{1)]^~^ symmetry to expose them. 
Accordingly, decompose the gauge field into diagonal or neutral(n) and 
off-diagonal or charged(c/i) parts A = A"' -\- A'^^. Gauge fix to a gauge such 
as maximal Abelian(MA) gauge [5] 

D;Af^d,Af-zg[A;,Af] = (4) 

which leaves a residual [t^(f)]^~^ symmetry. In MA gauge the residual 
[[/(I)]^""^ symmetry is 

/exp-^ \ ^ 

f^residual = '-. 5 2^^i=0. (5) 

V exp''^^ / '=^ 

Under this symmetry the N neutral fields {A'^)ii transform as photon fields 
whereas the N{N — I) charged fields {A'^^)ij transform as charged matter 
fields: 

(Alh ^ {A;)u + ^d,u,, {Al%, ^ {Al%, exp'-(---^) z, j G [I, N]. (6) 

Each {A'^^)ij carries two different U{1) charges so that the A'^^ fields induce 
interactions between the N photon species. Since Z]j=i(A")ji is an invariant 
under (6), an irreducible representation of [[^(f)] 
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oi^iA;h-j.j:iA;h. (7) 

8 = 1 

These angles transform under [t^(f )]^~^ as ^1 — ^ ^1 + -d^uJi and obey con- 
straint 

N 

E^; = o. (8) 

8 = 1 



The 6^ fields are compact angular variables since they are constructed 
from matrix elements of an SU{N) gauge field. As such, they may poten- 
tially contain monopole currents just like the compact U{1) field of CQED. 
If 't Hooft's hypothesis is correct, these would-be monopoles are MA gauge 
manifestations of nonAbelian gauge field features responsible for QCD con- 
finement. 

This procedure where QCD monopoles are identified using only the di- 
agonal 6^ components of the nonAbelian gauge fields is known as Abelian 
projection [6]. The quantum dynamics of the N 6^ fields is known as Abelian 
projection QCD or APQCD, whose action will be denoted Sapqcd- For- 
mally, Sapqcd is obtained by integrating out the A^^ fields in MA gauge 
from the QCD action: [7] 

Sapqcd ^ - log{ f[dA^'] eMSgcD) S[D;Af]}. (9) 

Sapqcd is a [t^(l)]^~^ invariant action depending on the N constrained com- 
pact 6^ fields. Numerically, Sapqcd is the action which would generate the 
importance sampling 6^ configurations currently made by Abelian projecting 
importance sampling QCD configurations. 

In this paper we will concentrate on the physical case, A^ = 3. A general 
U{1) X U{1) action consistent with APQCD symmetries is^ 

oo 3 3 

- Sapqcd = j:j:{j:nQpiL))+ E ^-(0^^, 0^p(l))} + " " " (10) 

L=l P(L) 8=1 «<J=1 

where P{L) refers to the square L X L plaquettes in the lattice and su- 
perscripts i and j refer to the three U{1) gauge fields. "• • •" refers to al- 
lowed operators we are neglecting, such as nonsquare Wilson loops, Polyakov 
loops(important at finite temperature), and nonlocal interactions between 



el). 



^The operators in (10) are not all independent. For example, by (8), cos 0p = cos(0p- 



Wilson loops of different shapes and sizes. We have also neglected explicit 
monopole operators, which are gauge invariant. Gauge, hermitian conjuga- 
tion of the SU{3) links{Ux, ij. —^ Ux,fM), and species permutation [8] symmetry 
require that 

• F{x + 27r) = F{x) and G{x + 27r, y) = G{x, y); 

• F[—x) = F[x) and G{—x^—y) = G[x^y); 

• G{y,x) = G{x,y). 

Therefore, up to a constant F and G must be of the form 

oo 

P{^) = J2l^i (^os{qx), (f f ) 

g = l 

OO OO 

G{x, y) = J2 J2{'^p,i cos(pa; + qy) + Bp^q cos{px - qy)j (f 2) 

p=l q = l 

where p, q are integers, Aq^p = Ap^q, and Bq^p = Bp^q. 

In (2) only Z = I, unit charge q = 1 monopoles k are elementary 
dynamical or "integration" variables because Scqed contains only L = 1 
plaquettes. L > 1 and q > 1 monopoles are composite(nonelementary) op- 
erators in cQED. Influenced by CQED, the most famous model of monopole 
confinement, previous studies of APQCD monopoles have been mostly re- 
stricted to L = q = 1. However, Sapqcd may in principle contain a whole 
range of nonzero fiq{L) (and also A^^g, Bp^q) coupling constants. Accord- 
ingly, extended L^ [9] and fractionally charged oc Ij q monopole currents [10] 
are potentially elementary dynamical variables in APQCD. The L^ and 1/ q 
monopole currents would arise from the BKT transformation of the corre- 
sponding cos qQp(L) operators in Sapqcd- Therefore, in this paper we extend 
the study oi L = q = 1 monopoles to q>l and L > 1 for a wide range of L 
and q combinations in APQCD. 



The possible relevance of extended and fractional monopoles is not just 
an abstract issue having to do with the form of Sapqcd] it also has a physical 
bearing. L > 1 monopoles have recently come of interest due to the open 
possibility that QCD monopoles are not pointlike and may have a nonzero 
physical core radius [11]. In this case it is reasonable to anticipate that the 
extended monopole operators of appropriate size L would couple better to 
the underlying physical objects than the usual L = 1 operator. Analogously, 
a similar argument can be made for q > 1 monopoles. There is no a priori 
way to know that q > 1 monopoles do not contribute in some quantitative 
way to the string tension or width of some representations of the QCD flux 
tube. 

Additionally, there is the issue of cross-species monopoles. We have so 
far been discussing "diagonal" monopoles, magnetic kinks in a single 6^ field. 
Diagonal monopoles correspond to nonzero /3g coupling constants. As defined 
in Section (3), "cross-species" monopoles are kinks in some linear combina- 
tion of two different 6^ fields. As illustrated in Ref. [10], such monopoles are 
elementary dynamical variables when corresponding Ap^q and Bp^q coefficients 
of G{x^y) are nonzero. Indeed, we will show that in APQCD cross-species 
monopoles have nontrivial behavior — in fact the same qualitative dependence 
on Z, q^ and lattice temperature r as the diagonal monopoles. Therefore, at 
least phenomenologically, in APQCD cross species monopoles are equally as 
interesting as their diagonal counterparts. 

Sections 2 and 3 present our numerical results. While we will provide 
some commentary, we do not claim to offer final conclusions or interpreta- 
tions. Section 2 discusses the correspondence between diagonal L > 1 and 
q > I monopoles and operators in Sapqcd- The nonAbelian gauge depen- 
dence of the Abelian projection and the temperature r dependence of the 
diagonal monopoles are described. We find that the diagonal monopole den- 



sity p in APQCD has its minimum ai L = q = 1. As explained, this result 
suggests fjq-ii^L = 1) is the dominant coupling in Sapqcd- We also find 
that the diagonal monopole current k has subtle but detectable temperature 
variations between and within the confining and finite temperature phases. 
These variations are such that the ratio of spatial to temporal monopole den- 
sities TZ is, as long as Z < l/r, an approximately Z-independent nontrivial 
indicator for the SU{'i) finite temperature deconfinement transition. For a 
range of L and q, TZ ^ 1 in the confining phase whereas TZ decreases very 
sensitively with temperature in the finite temperature phase. The redundant 
behavior of TZ at different L and q can have two possible origins: either the 
higher L and q monopoles are kinematically correlated to the L = q = 1 
monopoles and there is no independent L > 1 and q > 1 dynamics going 
on, or the L > 1 and q > 1 currents are acting independently and they only 
appear to be correlated because they share common dynamical rules. At 
this time it is not clear which of these explanations apply; we suspect it is a 
combination of both. 

Section 3 presents some exploratory results for cross-species monopoles 
as a function oi L, q, and r. For a range of Z, the ratios of spatial to tempo- 
ral cross-species monopole densities are nontrivial indicators for the SU{'i) 
deconfinement transition. Also, they vary sensitively with temperature in 
the finite temperature phase. 

2 Diagonal Monopoles 

In Toussaint-Degrand notation [12], define the generalized monopole current 

k,{L,q) = ^ E {(?0p(L))mod27r}, 9 = I,2,--- (13) 

^^ p(L)eC(L,^l) 



where C{L^ fi) refers to the L^ cube oriented in direction fi assuming D = 3+1 
dimensions. When q = 1 and L > 1^ k reduces to the Type I extended 
monopoles of Ref. [9]. Integer current k is topologically conserved for all 
L and q provided one uses an extended derivative when L > 1. Because 
k4 carries magnetic charge > l/qe if the fundamental representation Wilson 
line carries electric charge Q = e/2, we say q > 1 monopoles are fraction- 
ally charged or "fractional." By Dirac's quantization condition, fractional 
monopoles of charge l/qe can only interact with electric charge Q > qe/2 
Wilson loops. 

Our monopole density is calculated as^ 

p{L,q) ^ {ML,q)\) = Y: \HL,q)\ / T. ^ ■ (14) 

X X 

configs configs 

As previously discussed, only L = q = 1 monopoles are elementary dynamical 
dynamical variables in CQED because Scqed^ given in Eq. (1) contains only 
L = q = 1 plaquettes. Figure lA depicts a plot of the monopole density in 
CQED as a function of q above and below the critical point. As shown, while 
pil^q > 1) monopoles are also suppressed, p is most greatly suppressed at 
L = q = 1. As q becomes greater than 1, p converges to 

PR = ^ = .46, (15) 

the monopole density when links are completely random [15]. Similarly, 
as L becomes larger the extended links making up k[L^ q) — which are su- 
perpositions oi L = 1 links — become more and more disordered. Hence 
p{L » l^q) ^ PR for all q. 



^The sum over x ranges over all dual lattice sites; the cubes of k are permitted to 
overlap. Thusly defined, '^Zx^p-i-^'l) vanishes identically on our periodic lattices. When 
L = 1 our normalization of p agrees with Ref. [13]. When L > 1 our normalization of p 
disagrees with Ref. [14] by a factor of L^ . Also, since our SU{'i) results presented below 
are fixed at /3 = 6.0 we do not need to convert p to physical units for interpretation. 
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Figure 1: lA depicts the (j'-dependence of the P monopole density in the two 
phases of CQED; IB, in a mixed U{1) model with action comprised oi q = 2 
and q = 5^ L = 1 plaquettes. The sohd horizontal lines refer to random-link 
monopole density pR = 7/15. The dotted and dashed lines are guides-to-eye. 
As illustrated, p{L^ q) equals pR unless it is disrupted by an operator in the 
action of related L and q. p is in lattice units. 
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In a model with L > 1 and/or q > 1 plaquettes in its action, monopoles 
of corresponding L and q become elementary dynamical variables. Figure IB 
depicts p as a function of q for the action 

- S^^^ed = XI /^2 COs(20p(i)) + /35 COs(50p(i)) (16) 

P(l) 

at jSs = .7 and for a range of /32- When /32 < 1-0, /j(1,2), /j(1,5) and to a 
lesser extent /j(1,3), /j(1,4), and /j(1,7) are suppressed relative to pR. When 
q >> 5^ p uniformly approaches pR. On the other hand, at /32 = 1-2 the model 
crosses over to a weak coupling frozen phase where p is greatly suppressed 
for a range of q. 

In the /32 ^ 1-0 phase, much of the (j'-dependence of p depicted in Figure 1 
can be understood in terms of strong coupling arguments. If the action is 
zero, the U{1) links are random and p = pR. If we turn on a representation 
q^ L X L plaquette in the action, then such a plaquette can dress k[L^ q) 
and drive p{L, q) below pR. If a second plaquette operator is turned on, 
then k can be dress by the two plaquettes separately and in combination. In 
Figure 1 /j(1,3), /j(1,4), and /j(1,7) are suppressed because representations 
(j' = 3, 4, 7 are products oi q = 2 and q = b. /j(1, 1) and /j(1, 6) are much less 
suppressed because q = 1 and q = & cannot be made from 2 and 5. As /32 
increases to 1.2, the strong coupling picture breaks down. 

The strong coupling picture predicts the following behavior for the L- 
dependence of p. If the action contains only plaquettes of size Lo then p{L < 
Lo^q) ~ PR since it is geometrically very difficult for larger plaquettes to 
dress smaller monopole operators. Similarly, p{L > Lo^q) would differ from 
PR only due to high order strong coupling graphs. Hence only p at Lo and 
qo associated with corresponding operators in the action (or simple products 
thereof) would be appreciably suppressed. 

All APQCD results presented in this paper are on 24"^ X T, /3 = 6.0 
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Figure 2: 2A and 2B depict the q and L dependence of the APQCD monopole 
density p in different gauges, p is given in dimensionless lattice units. In the 
three gauges examined, the q = L = 1 monopoles are disrupted from the 
random pR = 7/15 value more than q ^ 1 oi L ^ 1 monopoles. Jackknife 
error bars are, in principle, drawn for all data points in this paper, although 
sometimes they are too small to be visible. 
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lattices with periodic boundary conditions. Temperature r is given in lattice 
units as T = l/T. On our lattices the N Abelian Polyakov loops {Pi) in 
APQCD vanishes when T = 40 and T = 8, and is nonzero when T = 6 and 
r = 4. Accordingly, the T = 40 and T = 8 lattices are below and T = 6 
and r = 4 are above the deconfinement temperature Tc. In the Figures, all 
quantities are plotted in lattice units with (sometimes unresolvable) jackknife 
error bars. The dotted lines are guide-to-eye lines; the solid horizontal line 
indicates the value of pR. 

In the confining phase of APQCD — corresponding to the confining phase 
of QCD [6, 16] — one can hope to obtain an indication of the dominant op- 
erators in S APQCD from p[L^q) in APQCD. Figures 2A and 2B are plots of 
/? at r = 40 in MA gauge and, for comparison. Landau and axial gauges. 
The Landau gauge over-fixes the gauge — it leaves no residual [t^(l)]^~^ sym- 
metry. This seems to "over-suppress" the monopole density relative to MA 
gauge for a reason we do not fully understand."^ 

In addition to what is shown, we have also evaluated p at other L and q 
combinations and, furthermore, produced similar results on 16"^ X 24, /3 = 5.7 
lattices. In all cases, we find p assumes its minimum value"* at L = q = 1^ and 
that p is a monotonically increasing function of L for fixed q and q for fixed 
L. Note that while Figure 2A is qualitatively similar to Figure lA, p{L^ q) 
in APQCD is quantitatively very different from its CQFD counterpart. For 
example, in MA gauge /j(1, 2)//j(1, 1) = 16. 0(. 0003) whereas this ratio in 
CQFD is 3.32(.004) at I3cqed = -99 and 15.3(.003) at /3cqed = 1-08. In 



^Note that one can always restore the [^7(1)]^""'^ symmetry by making Landau-gauge- 
violating random [?7(1)]^~^ gauge transformations. This yields an Abelian projection 
with the same monopole currents as in Landau gauge because the Landau gauge monopole 
currents are invariant by construction under [^7(1)]^""'^ gauge transformations. 

■*If we let integer q take on fractional values, the minimum nonzero value of p occurs at 
g = i corresponding to magnetic charges ±2. Since lattice monopole charges don't exceed 
magnitude 2, p{L, g < i) = 0. 
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agreement with an earlier analysis based on plaquette spectral densities [7], 
this suggests that between /3 = 5.7 — 6.0 the biggest coupling in Sapqcd is 
l3g-i(^L = 1). The last statement has two caveats. Firstly, we certainly do 
not rule out the existence or significance of other fiq{L) coupling constants. 
Secondly, any substantial nonzero Ap^q or Bp^q coefficients in Sapqcd will 
infiuence p[L^q)^ possibly obscuring our interpretation. 

Figures 3 and 4 depict the temperature r dependence of q = 1 monopoles 
in MA gauge for a range of L. q > 1 monopoles (not shown) behave simi- 
larly although, since p{L^q > 1) is closer than p[L^l) to pR^ they are less 
sensitive indicators than q = 1 monopoles. While Figure 3 may seem to 
merely extend the L = q = 1 SU{2) results of [f 7] to Z > f in SU{3)^ let us 
mention the following technical distinction which is nontrivial. In Ref. [17] 
the high temperature phase was achieved by varying fiqcD on a lattice with 
a fixed number of timeslices. Varying fiqcD has the disadvantage that it 
simultaneously changes the physical size of the I"^ monopole operator. Thus, 
this method introduces a hidden variable which may (or may not) have an 
effect on the observed monopole density. On the other hand, in this paper 
we fix fiqcD = 6.0 and achieve finite temperatures by varying the number 
of timeslices T . Since our physical lattice spacing is fixed, the physical size 
of our I"^ monopole operator is fixed and we are offering a truly independent 
verification of Ref. [17] 's original result. 

Figures 3A and 3B depict the behavior of the spatial and temporal 
monopole densities for Z = I — 4. The temporal density is p of Fq. (14); 
the spatial densities are analogously defined using the spatial A;^ compo- 
nents. The temporal and spatial densities are only mildly r-dependent with 
the L = 1 densities being the most sensitive. As was noted in Ref. [17] 
for L = q = 1 SU{2) monopoles, contrary to naive expectations APQCD 
monopole densities do not decrease dramatically across the deconfinement 
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Figure 3: 3A and 3B show that the spatial and temporal monopole densities 
both decrease across the deconfinement transition and both increase within 
the confined phase with increasing r. Within the deconfined phase the spa- 
tial density decreases whereas the temporal density rises with increasing r. 
Therefore, in 3C the spatial-temporal asymmetry ratio TZ has enhanced sen- 
sitivity to T in the deconfined phase for all L. The 4"^ monopoles suffers 

boundary effects at r = f /4. 
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Figure 4: 4A and 4B present the q and Z dependence of the ratio TZ. While 
TZ converges to the uninteresting value of unity at q > A at all temperatures, 
it is only very mildly Z-dependent. At this time we do not understand 
the significance of such near Z-independence, which indicates that certain 
features of Z > 1 monopoles correlate to temperature as positively as Z = 1 
ones. 
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transition as occurs in CQED. While both densities rise with r within the 
confinement phase, the spatial density falls while the temporal density mildly 
rises with r in the high temperature phase. 

There are two comments to be made here. Firstly, monopoles are needed 
at finite temperature to account for the string tension of spatially oriented 
Wilson loops [18]. Hence, consistency with their confinement-causing role 
requires that the monopole currents do not disappear at finite temperatures. 
Secondly, the Abelian Polyakov loop Pi = exp i j 0\ gains a nonzero expec- 
tation value at finite temperatures because the time-like link angles freeze 
(modulo U{1) gauge transformations) to 6\ — )• 0. Since the spatial compo- 
nents of the monopole currents share time-oriented links with Pj-, and small 
link angles tend to decrease the probability of forming a monopole kink, the 
falloff of spatial densities and the rise of {Pi) probably share a common ori- 
gin: freezing of the time-oriented links with increasing r. Note that k4 does 
not contain time-oriented links. 

Figure 3C depicts the ratio TZ of spatial to temporal monopole densities. 
As shown, 7^ is a sensitive finite temperature indicator not just at L = 1 [6] 
but for all L = 1 — 4:. 

Figures 4A and 4B reveal the q and L dependence of TZ at different 
temperatures. While TZ becomes less and less informative as q increases, 
what is compelling is that 7Z is approximately L- invariant as long as L is 
smaller than the smallest lattice width T. 

3 Cross Species Monopoles 
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Figure 5: While cross-species monopoles are generally denser than the diag- 
onal ones, they have the same general q^ L^ and temperature dependence as 
the latter. As depicted in 5A the cross-species monopole density is smallest 
at Z = 1 and rises quickly to pR at bigger L for all r. 5B shows that the 
asymmetry ratio 7^'^™'^'^ is sensitive to temperature in the finite temperature 
phase. Note the boundary effect on the L = A current at r = 1/4. 
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In addition to k[L^ q) one can also consider cross-species monopole currents 
defined by 

k--{L,p,q,±)^^ E {(p0}.(L)±?0^(L))mod27r} (17) 

^^ P(L)eC(L,^) 

for integers p and q. Cross species monopoles are elementary dynamical 
variables arising from the interspecies interaction operator G of Eq. (12). 
For a special case of (10) we illustrated in Ref. [10] how both k and k'^^'^^^ can 
simultaneously occur as elementary dynamical variables. In such U{1) X U{1) 
models the dual Meissner effect depends on the combined status of k and 
^cross^ which may in principle condense or freeze out independently for each 
value of L and q. 

In APQCD k'^^°^^[L^q^q^-\-) is equivalent to k[L^q) by constraint (8). 
Henceforth we focus on k'^^°^^[L^q^q^—)^ for which some MA gauge results 
are presented in Figure 5. p'^^"^^ and 7^'^™'^'^ are defined analogously to their 
diagonal counterparts. Our exploratory calculations indicate that generally 
^cross^ while denser, has qualitative features reminiscent of k. As demon- 
strated in Figure 5B, 7^'^™'^'^ is also a nontrivial indicator for the SU{3) finite 
temperature transition. However, it is noticeably more Z-dependent than its 
diagonal counterpart. 

We stress that it is unclear at present whether in APQCD k'^^"^^ is truly 
dynamically independent of k. 
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